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METHOD AND APPARATUS FOR AUTOMATIC DETECTION AND 
IDENTIFICATION OF A BROADCAST AUDIO OR VIDEO 
PROGRAMMING SIGNAL. 

5 

BACKGROUND AND SUMMARY OF THE INVENTION 

This invention relates to the automatic detection and identification of broadcast 
programming, for example music or speech that is broadcast over radio, television or 

10 the Internet, or television signals, whether broadcast as analog, digital or digital over 
the Internet. By "Broadcast" it is meant any readily available source of content, 
whether now known or hereafter devised, including, for example, streaming, peer to 
peer delivery of of downloads or streaming or detection of network traffic comprising 
such content delivery activity. The system initially registers a known program by 

15 digitally sampling the program and separating the digital sample stream into a large 
set of short segments in time. These segments are then processed to extract particular 
feature sets that are characteristic of the segment. The invention processes each set of 
features to produce a numerical code that represents the feature set for a particular 
segment of the known program. These codes and the registration data identifying the 

20 program populate a database as part of the system. Once registration of one or more 
programs is complete, the system can then detect and identify the presence of the 
registered programming in a broadcast signal by extracting a feature set from the 
input signal, producing a numerical code for each time segment input into the system 
and then comparing the sequence of detected numerical codes against the numerical 

25 codes stored in the database. Various testing criteria are applied during the 
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comparison process in order to reduce the rate of false positives, false negatives and 
increase correct detections of the registered programming. The invention also 
encompasses certain improvements and optimizations in the comparison process so 
that it executes in a relatively short period of time. 

BRIEF DESCRIPTION OF THE DRAWINGS 
The components of the media broadcast monitoring system. 
An illustration of the data flow of the detection algorithm from a series 
of frames of an audio program to detection of the program's identity. 
The flowchart of the Pattern Generation Module. 

Example of how original frequency band boundaries lead to pattern 
mismatches between the original frame signatures and the signatures, of 
the same audio program played at a faster speed. 

Example of how changing the frequency band boundaries yields an 
improved match between frame signatures of the original audio 
program and the same audio program played back at fast and slow 
speeds. 

The new frequency band boundary setting leads to robustness of the 
audio detection algorithm even with +/-2% speed variations in the 
audio program. 

The schematic of the DBS operation flow. 
The flowchart of the SRR Algorithm. 
Example Calculation of Frequency Band Boundaries 



Figure 5: 



Figure 6: 



20 Figure 7: 
Figure 8: 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Background. 

5 The present invention relates to the automatic recognition of widely disseminated 
programming, such as radio, television or digitally delivered content over the Internet. 

Owners of copyrights in broadcast programming, including advertisers, need to 
measure when and where their programming has been broadcast in order to correctly 

10 compute performance royalties, confirm compliance with territorial restrictions or 
verify that certain advertising has been aired as scheduled. The traditional method for 
monitoring the radio or television has involved using humans to listen or watch and 
then record that which they hear or see, or alternatively, rely on the broadcast records 
of radio and television stations. This is a labor intensive process that has limited 

15 efficiency or accuracy. It is an object of the invention to use advanced computing 
systems to fully automate this process. In this manner, audio or video content is 
registered into the system, and then, in the case of audio detection, radio, the 
soundtrack from television or other sources of widely distributed audio content are 
input into the system. In the case of video, the video signal is input into the system 

20 from whatever its source. By means of the invention, the detection and identification 
of registered programming content takes place automatically. 

Prior Art: 



WO 2005/081829 PCT/US2005/004802 

-4- 

A number of methods have been developed to automate the detection of broadcast 
programming. These techniques generally fall into one of two categories: cue 
detection or pattern recognition. The cue detection method is exemplified by U.S. Pat. 
Nos. 4,225,967 to Miwa et. al.; 3,845,391 to Crosby and 4,547,804 to Greenberg. 
5 These techniques rely on embedded cues inserted into the program prior to 
distribution. These approaches have not been favored in the field. In audio, the 
placement of cue signals in the program have limited the acceptance of this approach 
because it requires the cooperation of the program owners and/or broadcasters«thus 
making it impractical. 

10 The pattern recognition method generally relies on the spectral characteristics of the 
content itself to produce a unique identifying code or signature. Thus, the technique 
of identifying content consists of two steps: the first being extracting a signature from 
a known piece of content for insertion into a database, and the second being extracting 
a signature from a detected piece of content and searching for a signature match in the 

1 5 database in order to identify the detected content. In this way, the preferred approach 
relies on characteristics of the broadcast content itself to create a signature unique to 
that content. For example, US Patent No. 4,739,398 to Thomas, etal. discloses a 
system that takes a known television program and creates for each video frame, a 
signature code out of both the audio and the video signal within that frame. More 

20 recently, similar detection systems have been proposed for Internet distributed 
content, for example application PCT WO 01/62004 A2, filed by Ikeyoze et. al. 

For audio by itself, U.S. Pat. No. 3,919,471 to Moon discloses an audio identification 
system where only audio signals are used, but it is of limited utility because it 
25 attempts to correlate an audio program represented by a limited time slice against the 
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incoming broadcast signal. The disclosed method of matching in Moon is highly 
compute intensive because it relies on direct signal correlation. Further, this approach 
is unfavorable because it has been found to be limited in accuracy, especially if the 
program is time compressed or altered in other ways prior to detection. It is also prone 
5 to false positive identifications and is computationally uneconomic if the size of the 
time slice is expanded to improve its correct identifications. Lert, et. al. describes in 
U.S. Pat. No. 4,230,990 a way to mitigate the computational workload of the 
correlation method by combining it with the coding method of the first category: 
either an artificial code or some other naturally occurring marker is detected in the 

10 program indicating the beginning of a section of the program, and then a feature 
signature is measured at a pre-determined amount of time later. This method has 
limited utility in audio-only applications, where either an audible code has to be 
inserted into the audio to create the cue, thus degrading it or requiring cooperation of 
the content source, or reliance on natural markers indicating the start of a new audio 

15 program which is highly unreliable. In U.S. Pat. No. 4,677,466 Lert, et. al. further 
describes an improvement on the invention that waits until a "stability condition" has 
occurred in the signal before measuring and calculating a signature, but the reliability 
of the method is limited by the size of the sample time slice. U.S. Pat. No. 4, 739,398 
to Thomas et. al. addresses the data processing load problem by randomly choosing 

20 portions of a signal to sample as input to the invention's signature generating process. 

U.S. Pat. Nos. 5,436,653 to Ellis, et. al. and 5,612,729 to Ellis, et. al., disclose a more 
complex way of calculating a unique signature, where the audio signature 
corresponding to a given video frame is derived by comparing the change in energy in 
25 each of a predetermined number of frequency bands between the given video frame 
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and the same measurement made in a prior video frame. However, the matching 
technique relies on a combination of the audio and video signatures or the use of a 
natural marker, in this case, the start or ending of a program. Thus, this method 
suffers the same problem as Lert with regard to audio-only programming. 

5 

In addition, U.S. Pat No. 5,918,223 to Blum, et. al., discloses the use of audible 
features within audio programming to create a single signature value for each audio 
program, particularly the group of amplitude, pitch (i.e. fundamental), bandwidth, 
bass (i.e. rhythm analysis), brightness (i.e. shape of the frequency response of the 

10 program), and Mel-frequency cepstral coefficients. The aggregation of these detailed 
features across long periods in the audio produce highly variable results, and do not 
possess sufficient robustness in real-world broadcast situations. U.S. Pat. No. 
5,210,820 and 4, 843, 562, both to Kenyon, discloses a digital circuit that uses the 
envelope (e.g loudness) features in the audio signal in order to create a signature. The 

15 approach is designed to address the time compression problem by application of time 
warping techniques. Reliance on loudness has other robustness problems that also 
make it difficult to use in real-world environments. U.S. Pat. Application No. 
20030086341 filed by Wells, Maxwell, et. al, discloses a system where an audio 
signature is created using pre-determined numbers of digital samples counted from 

20 pre-determined locations from the start point of the music. This approach is much less 
reliable for broadcast or cases where the audio is detected in analog form, or in cases 
where the playback of the programming has changed speed, frequency equalization 
from the original track has been applied, or the audio dubbed into the programming 
segment. 
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The present invention describes a system and method whereby identification of 
known audio or video programming can be done without any reliance on a tandem 
video signal (in the audio case) or normative markers in the signal indicating a known 
time in the program and with unique and novel ways to calculate codes representing 
5 the characteristics of the audio program without requiring impractical computational 
capabilities. Benefits of this system and method are accuracy, speed, robustness to 
playback speed variation and the ability to perform the identification process in real 
time, without reliance on any embedded cue or watermark. In addition, the present 
invention takes advantage of the availability of low cost, high performance computing 
10 platforms in order to implement a high speed database searching methodology. 



Detailed Description. 
A. Overview 

The broadcast monitoring and detection system embodying the invention works in 
15 two phases: registration and detection. During the registration phase, known 
programming content is registered with the system by sending the program, as digital 
data, into the system. A series of signatures, in the case here, a pattern vector and 
more generally in the art a "fingerprint" or "signature", are stored as a sequence of 
data records in a database, with the identity of the program content cross-referenced to 
20 them as a group. During the second phase, unidentified programming is input into the 
system. Such programming can include radio, television, internet broadcasts or any 
other source of audio or video programming, whether terrestrial broadcast, satellite, 
internet, cable television or any other medium of delivery, whether now known or 
devised in the future. While such programming is being monitored, the pattern 
25 vectors of the programming (or any other signature generating technique) are 
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continually calculated. The calculated pattern vectors are then used to search for a 
match in the database. When a match is found and confirmed, the system uses the 
cross-referenced identity in the database to provide the identity of the content that is 
currently being played. In the preferred embodiment, the system is software running 
5 on a computer, however, it is envisioned that special purpose hardware components 
may replace parts or all of each module in order to increase performance and capacity 
of the system. 



In the preferred embodiment, a computer containing a central processing unit is 

10 connected to a sound card or interface device into which audio programming is 
presented. During the registration phase, the CPU fetches the audio or video data 
from the sound card, calculates the pattern vector data, and then, along with timing 
data and the identity of the program, these results are stored in a database, as further 
described below. Alternatively, the data may be loaded directly from authentic 

15 material, such as compact discs, mp3 files or any other source of digital data 
embodying the signal. For non-audio applications, the source of material can be DVD 
disks, masters provided by movie studios, tapes or any other medium of expression on 
which the program is fixed or stored. Of course, for some material which may not 
have a readily available source, then the audio or other program signal is used in the 

20 following manner. If the system periodically detects an unknown program but with 
the substantially the same set of signatures each time, it assigns an aribitrary identifier 
for the program material and enters the data into the database as if the program had 
been introduced during the registration phase. Once the program identity is 
determined in the future, then the database can be updated to include the appropriate 

25 information as with authentic information while at the same time providing the owner 
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of the programming the use data detected even when the identity of the program was 
not yet known. The database, which is typically a data file stored on a hard drive 
connected to the central processing unit of the computer by means of any kind of 
computer bus or data transmission interface, including SCSI. 
5 During the detection phase, the CPU fetches the program data from the sound card or 
video card, or loads it from a data file that may be stored on the computer hard drive 
or external media reader. The CPU calculates the pattern vector data, and then, along 
with the timing data, submits database queries to the database stored on the hard 
drive The database may be the same hard drive as in the computer, or an external 

10 hard drive accessed over a digital computer network. When matching data is found, 
the CPU continues to process the data to confirm the identification of the 
programming, as described further below. The CPU can then communicate over any 
of a wide variety of computer networking systems well known in the art to deliver the 
identification result to a remote location to be displayed on a screen using a graphical 

15 user interface, or to be logged in another data file stored on the hard drive. The 
program that executes the method may be stored on any kind of computer readable 
media, for example, a hard drive, CD-ROM, EEPROM or floppy and loaded into 
computer memory at run-time. In the case of video, the signal can be acquired using 
an analog to digital video converter card, or the digital video data can be directly 

20 detected from digital video sources, for example, the Internet or digital television 
broadcast. 

The system consists of four components. Figure 1 shows the interconnection of the 
four modules: (1) a signal processing stage at the front end, (2) a pattern generation 
module in the middle, (3) followed by a database search engine module, and (4) a 
25 program recognition module at the end. During the registration phase, the results of 
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the pattern generation module, which creates signatures for known audio or video 
content, are stored in the database and the search and pattern recognition modules are 
not used. 

5 The function of each module is described in further detail below: 
1 . Sound Acquisition (SA) Module 

The SA module, (1), receives audio data from a sound detection circuit and 
makes it available to the remaining modules. Practitioners of ordinary skill 

10 will recognize that there are a variety of products that receive analog audio or 
video and convert those signals into digital data. These devices can be any 
source of digital audio data, including an interface card in a personal computer 
that converts analog audio into digital audio data accessible by the computer's 
CPU, a stand alone device that outputs digital audio data in a standard format 

15 or a digital radio receiver with audio output. Alternatively, pre-detected signal 
in digital form can be accessed from storage devices connected to the system 
over typical data networks. The SA module regularly reads the data from the 
digital interface device or data storage and stores the data into a data buffer or 
memory to be accessed by the Pattern Generation module. Practitioners of 

20 ordinary skill will recognize that the typical digital audio system will provide a 
digital word at regular intervals, called the sampling rate. The sequence of 
digital words representing the audio signal are the digital audio samples. The 
invention organizes the samples into a series of time frames, which consist of 
a predetermined number of samples. The time frames are stored in sequence. 

25 Alternatively, data structures, stored in the computer memory (which includes 
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the hard drive if the operating system supports paging and swapping), may be 
used where the time frames are not physically stored in sequence, but logically 
may be referenced or indexed in the sequence that they were detected by 
means of memory addressing. 

5 

In the preferred embodiment, the audio signal is conditioned in a manner 
known in the art, including low-pass filtering. In the preferred embodiment, 
the signal is sampled at a rate of 8000Hz within the SA Module. In the 
preferred embodiment, 16,384 samples constitute a single frame. At this rate, 
10 the signal must be lowpass filtered for anti-aliasing purpose before being 
sampled. Higher sampling rates may be used, however, with the appropriate 
adjustments in the downstream calculations, as explained below. 

In the case of video programming, the sound acquisition module essentially 
acts in an analogous manner: the video signal is acquired as a digital video 

15 signal, and converted to the frequency domain using well known methods on a 
video frame by frame basis. The invention will be explained in detail as 
applied to audio through a description of the preferred embodiment. However, 
the system and processes described are applicable to video as well as audio, 
where a signature or pattern vector has been periodically derived from the 

20 video signal. Reference is made to "A Technical Introduction to Digital 
Video", Charles A. Poynton, John Wiley & Sons, New York, © 1996. 
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2. Pattern Vector Generation (PG) Module 

The PG module operating during the detection phase, (2), fetches the stored 
digital audio or video samples that were detected and stored by the SA 
Module. Once a frame of the samples is received, the PG module will 
5 compute the pattern vector of the frame and, when in detection phase, send the 
pattern vector to the Database Search Module in the form of a database query. 
During the registration phase, the PG module calculates the pattern vector in 
order that it be stored in the database, in correlation with the other relevant 
information about the known audio or video program. The calculation of the 
1 0 pattern vector is described further below. 



Inter-Frame Distance. 

For each incremental audio sample, a new frame can be started. That is, each audio 
sample may be the constituent of N overlapping frames when N is the number of 

15 samples in a frame. The distance between these overlapping frames is the inter-frame 
distance. The shorter inter-frame distance for pattern generation mitigates the problem 
of program start-time uncertainty. Shorter inter-frame distances produce better 
results when the start time is unknown. In the preferred embodiment, the value of 
4,000, around V4 of a frame, is used during the audio program registration phase. 

20 Other distances may be used either to increase accuracy or reduce compute time and 
storage overhead. Thus, in the preferred embodiment, the first frame in the database 
of known audio programs corresponds to audio samples 1 to 16,384, the second 
corresponds to samples 4001 to 20,384, and so on. During the detection phase, the 
inter-frame distance is set to be equal to one frame length. Thus, the first frame of the 
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detected audio program contains samples 1 to 16,384, the second frame contains 
samples 16,385 to 32,768, and so on. 

Even though the uses a preferred embodiment setting of sampling rate of 8000Hz, 
frame-size of 16384 samples, inter-frame distance of 4000, a different sampling rate 

5 may be used with varying results. For example, for a sampling rate of 16000Hz 
(double the preferred setting), results in a frame number size of 32768 (double in size 
but the same in time duration), inter-frame distance of 8000 (inter-frame distance is 
the same at 0.5 sec) and generates almost identical pattern vectors as when using the 
preferred settings. The only further change is to determine which Fourier Transform 

10 (FFT) coefficients would be included in each sub band used to calculate the pattern 
vectors. For example, with the preferred settings, (ignoring the speed compensation 
scheme explained below), band 1 comprises the 66th to 92nd FFT coefficients. Then 
with the alternate example above, the FFT coefficients will be the 32nd to 94th. The 
calculation of the pattern vectors, which is presented assuming the sampling rate of 

15 8000 Hz, is adjusted accordingly. 

In the case of video, the pattern vectors are derived from the two-dimensional FFT 
transform of each frame of video. The video frames can be considered analogous to 
the samples in audio. Thus the vertical and horizontal FFT coefficients can be 
20 collected across the video frames to build pattern vectors for each time frame, the 
time frames constituting a group of video frames. Practitioners of ordinary skill will 
recognize that the approaches may be combined, in that features of the audio 
soundtrack of the television program can be combined with features of the video 
signal of the same program to produce the pattern vectors. 
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3. Database Search (DBS) Module 

Upon the reception of a query generated by the PG module, this module, (3), will 
search the database containing the sequence of pattern vectors of known programming. 
If a match is found, then the module returns a set of registration numbers otherwise 
referred to herein as program-id's and frame-id's, referred to also as frame numbers, 
corresponding to the identities of a set of audio or video programs and the time frame 
numbers within these programs where the match occurred. If the search of the database 
fails to find a match, the DBS Module will issue a NO-MATCHED flag. It is 
contemplated that aspects of the invention for the DBS Module are applicable to any 
kind of data set containing signal signatures, even signatures derived using techniques 
distinct from those used in the Pattern Vector Genreation module. 

4. Program Detection and Identification (SDI) Module 

This module, (4), constantly monitors the matching results from the DBS on the most 
15 recent contiguous of N time frames, as further described below. In the preferred 
embodiment, N is set to five, although a larger or smaller number may be used with 
varying results. Two schemes are used to determine if any audio or video program 
has been positively detected. The first is a majority voting scheme which determines 
if, within each thread of matching pattern vectors among N, the number of frames that 
20 possess a valid sequence pass a designated majority of the block of frames. The 
second is a frame sequencing scheme which follows each of the potential thread and 
counts how many frames within that thread constitute a valid sequence. If there exists 
a thread(s) where a majority of the contiguous frames satisfy the frame sequencing 
requirement, then the program (whether audio or video) is deemed detected in that 



5 



10 
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thread. Either or both schemes are used to suppress false positive detections and to 
increase the correct detections. In the preferred embodiment, both schemes are used. 
Given a program (or more than one) that is detected, the SDI module will initiate two 
modes: 1. Identification mode : in this mode, the module logs all the reference 
5 information of the detected program, including title, songwriter, artist, record label, 
publishing company or any other information input during the registration phase of 
the system, along with the time when the program is detected, and the time into the 
program that the detection was made. This information will be registered on the 
detection log. 2. Tracking mode : In this mode, the module tracks each detected 
10 program by monitoring if the queried result of every new frame of the broadcast is 
obeying the sequencing requirement, described below. The algorithm is locked in this 
mode until the queried results cannot be matched with the sequencing requirement. 
Upon the exiting from the tracking mode, a number of detection attributes, including 
the entire duration of the tracking, and the tracking score, will be logged. 

15 

The pattern vector generated by the PG Module is sent to the DBS Module in order to 
conduct a search of the database for a match. The output is either a NO-MATCHED 
flag, which indicates that the DBS fails to locate a frame within the database that 
passes the search criteria; or the program-id's and frame-id's of the library patterns 

20 that pass the search criteria. 

The SDI Module collects the output from the DBS Module to detect if a new 
audio program is present. If so, the detected song is identified. Figure 1 is an 
illustration of the flow of the algorithm from a frame of audio to its result after 
detection. With regard to the application of the invention to video, the 

25 operation is analogous, once the pattern vectors have been generated. It is 
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contemplated that aspects of the invention for the SDI Module are applicable 
to any kind of data set containing signal signatures, even signatures derived 
using techniques distinct from those used in the Pattern Vector Genreation 
module. 

5 

Pattern Vector Generation. 

The PG module reads in a frame of signal, preferably consisting of 16,384 samples, 
with sampling rate preferably set at 8,000 samples per second. Thus, the frame length 

10 is approximately two seconds in time. More or less samples or frame widths in time 
may be used with varying results. Given x = [x x x 2 x 16384 ], the vector 
containing a frame of signal, where each x f is the value of the nth audio sample, an N 
element pattern vector is calculated with the following steps. In the preferred 
embodiment, N is equal to 31. Practitioners of ordinary skill will recognize that the 

15 value of N is arbitrary, and can be increased or decreased with varying results. For 
example, decreasing N reduces the compute time and memory requirements, but may 
reduce accuracy. Increasing N may do the opposite. Also, the method presented will 
assume that a 31 element pattern vector is being used in the calculation in order to 
simplify the presentation of the invention. Practitioners of ordinary skill will 

20 recognize that the same methodology will work when N is increased or decreased, 
depending on whether the goal is increased accuracy or reduced computer complexity. 



1. The Fourier transform of x is calculated with the number of points equal 
to the number of samples in the frame, in order to get the spectrum vector 
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SOQOsamples I sec A . ool7 
The spectral resolution of the transform is = ^ 384sample — ^" 



ies 



Segregate the FFT spectral values into frequency bands of a specified width, where in 
5 the preferred embodiment, the width is 64 Hz. The invention will be further 
explained in terms of the preferred embodiment in order to simplify the presentation, 
but without limitation to the extent of the invention claimed. 

Band #1 is from 0 to 64Hz, Band #1 encompasses FFT coefficients Xi 

tO Xl31 

10 Band #2 is from 64 to 128Hz, Band #2 encompasses X 132 to X 26 2, and 

so on. 



2. Compute the centroid (or center-of-gravity COG) of each band: 
15 m 



20 



Pk 131 



In the preferred embodiment, only Band 2 to 32 is used because Band 1 is 
the lowest band including zero Hz, which is normally not useful in FM 
radio transmission; and Band 32 covers the band up to 1,800Hz, which is 
typically sufficient bandwidth to encode a fingerprint of the audio. Of 
25 course, higher bands or lower bands can be used if required. The inclusion 

of higher or lower bands to account for signal characteristics can be 
determined empirically. The first step, where the FFT coefficients are 
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collected in order to calculate the centroid in step 2 is different in the case 
of video. Ih the video case, the FFT coefficients have to be selected from 
locations either in the complex plane or on the 2-dimensional spatial 
frequency plane as described on page 23 of Poynton, incorporated herein 
5 by reference. These locations are analogous to the frequency bands on the 

audio case. In a manner analogous to using predetermine frequency bands 
in audio, predetermined regions on the vertical/horizontal plane in the 
frequency domain can be defined and the FFT coefficient values in each 
regions used to calculate an element corresponding to that region. Once 
10 this selection is made, the centroid can be calculated in an equivalent 

manner. It is advantageous to ignore the frequency region encompassing 
the frame rate, sync rate, subcarrier, or line rate. The end result is 
essentially equivalent to the case of audio: that each time frame of video 
will have a pattern vector associated with it that is stored in a database. 

15 

After Step 3, a 31-element vector is obtained: 
c = [p 2 p 3 p 32 \= z [ c i c i c 3i]- hx the preferred embodiment, a further 

step converts c to an unsigned integer. The unsigned format is used because all the 
elements in c are positive in the interval of (1, 131). A further calculation on c 
20 normalizes each element to a value between 0 and 1 by exercising the division by 
131, the number of FFT components within each band: 

0<c, =-^<l 
' 131 

In the preferred embodiment, each element is then converted to the unsigned 16-bit 
integer format for convenient storage and further processing. In order to decrease the 
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compute time downstream, each FFT coefficient or c f is tested relative to a minimum 
threshold value. The downstream processes are set to ignore these elements, for 
example, by not including these elements in downstream sets that are collected for 
further calculation. Figure 3 shows a flowchart of this module. In the preferred 
5 embodiment, both the FFT in step 1 and the centroid (COG) computation in step 3 are 
typically calculated using double precision floating point instructions. 



Speed Compensation Scheme 

Practitioners of ordinary skill in the art will recognize that for a variety of reasons, 
10 broadcast programming is often sped up from the speed of the original programming. 

Therefore, it is critical that any automatic audio program detection system be robust 

when the detected audio program may differ from the speed of the audio provided 

during the registration phase. In order to alleviate this problem, a modification to the 

pattern vector generating formula is used: 
1 5 (a) The modification is to have a different number of FFT components (i.e. 
bandwidth) of each band in step 2. 

(b) In the preferred embodiment, the modification to the pattern vector generation 
formula is only applied to the incoming broadcast audio signal during the 
detection phase, not to the pattern generation process applied during the 
20 registration phase of the audio program. Practitioners of ordinary skill will 

recognize that the use of the alternative frequency bands described above for 
the detection phase can alternately be performed during the registration phase 
with substantially the same result. 

The specific detail of this modification is described below: 
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The formulation is based on the scaling property of the Fourier Transform. 
A time speed up version of a song is a time-scaled version of the original: 
x(t) — speedup > x{at) ; a > 1 where a is the rate of speedup and x(t) is the detected 
sample at time t Note that for a > 1, the time axis is "compressed". If the song is 
5 sped up by 2%, we have a = 1.02 . 

With the scaling property, the factor a can be used to adjust the values of the Fourier 
Transform: 



Thus, the spectrum of a fast playback, or speedup version of a song is stretched. With 
10 a 2% speedup rate, the Fourier Transform frequency component at 100Hz without any 
song speedup, is shifted to 102Hz after speedup. This implies that, if there exists a 
2% speedup rate in the detected song, the bandwidth in step 2 should be adjusted 
accordingly to 1.02 x 64Hz = 65.28Hz, and hence the number of FFT components 
within each band should be adjusted to the roundoff of 131 x 1.02, which is equal to 
15 134. There are two formulae to calculate the amount of FFT components in each 
band, both based on the original number of FFT components, which is equal to 13 1. 

Formula 

(1) Given the speedup rate r. 
20 Start at Band #1, which encompasses FFT coefficients Xi to X z( i) , where z(l) = 



roundoff of 13 lx (1+r). 

(2) Compute iteratively each z(k) = roundoff of [z(k - 1) + 131x (14r)] for k = 2 to 
32. Band #m consists of FFT coefficients of X z (m-i)+i to X z(m ) . 




■+*(/) 
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(3) Compute the centroids (COG) from Band #2 to #32 with the new band 
partitions calculated above. Exercise the normalization by dividing each centroid 
(COG) by the number of FFT components in the corresponding band. 

5 The difference with and without the compensation is shown in Figures 4 and 5. 
Figure 4 shows Original band setting leads to pattern mismatches between the original 
and its speedup variant. Figure 5 shows that the modified band setting yields very 
good pattern matching behavior given that the speedup rate is known. 

1 0 Robust Pattern Vector Generation Formula 

The pattern vector generation formula described above can be further refined in order 
to provide robust matching. This refinement may also be used instead of the prior 
formulation. Besides causing the frequency axis to stretch, another effect of speedup 

15 is the shift of the boundaries in frequency of every band. The refinement is to 
compensate the shift of the boundaries of a band by extending the width of the band, 
such that the amount of the shift due to playback speed is within a small percentage 
compared with the band width. Thus, there is no modification of the algorithm - that 
is, calculating centroids as pattern vectors - except that the band locations are 

20 changed. The modified band boundaries are used during the registration process to 
create the stored pattern vectors. Practitioners of ordinary skill will recognize that 
several alternative methods may be used to calculate frequency band widths that 
exhibit the same property, that is, extending the band width such that the frequency 
shift due to playback speed variation is comparatively small, where the percentage 

25 frequency shift due to playback speed changes is a small percentage of each 
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frequency band width. Further, it is contemplated that this technique will work for any 
method of calculating a signature in the signal that is based on segregating the FFT 
coefficients into frequency bands. One method to calculate modified band boundaries 
that exhibit this effect is described below as the preferred embodiment. 

5 

Algorithm to compute new band boundary locations: 

Let the starting and ending indexes of band number k in the frequency domain be s k l 
and s k2 respectively, that is the index of the FFT coefficients. For example, index 
s u is equal to 1 and corresponds to the first FFT coefficient for 0 Hz. A shift-to- 
10 bandwidth ratio is assumed, which is the expected maximum speedup in percent 
divided by the percentage of bandwidth that the shift should not exceed. In the 

preferred^ 

in order to increase accuracy or reduce compute complexity. 

1. Start from band k=l, whose starting location s ui = 1 . Assuming a 2% 

1 5 speedup, the location is shifted by 0.02 to 1 .02, which after roundoff is still 

equal to L Roundoff is necessary because the result indices must be integers. . 
Assuming the shift-to-bandwidth ratio to be equal to 0.4 (which is 2% shift 
divided by 5% bandwidth, the amount that shift should represent) of the 
bandwidth of Band #1, then the ending location s x 2 = (1 + .02/. 05) x s u =1.4, 

20 or 1 after round-off. 

2. Now proceed to compute the two locations for Band #2. The starting location 
s 2l =2. Given 2% shift and 5% shift-to-bandwidth ratio, we obtain s 2 j2 =3. 

3. Continue the iteration until all the FFT components are exhausted. In the 
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preferred embodiment, the result (both lower order bands, , s k l < 64 , 
corresponding to 3 1 .25 Hz, and higher order bands, s k l > 5500 , 

corresponding to 2,686 Hz, are not used. 

4. When k equals 9, then s 9 2 = 66 , and when k= 10, s\q 9 \ = 67 and so on. 

5 In order to avoid overflow because the bandwidth of each band along k increases 

exponentially with k, the preferred embodiment has set arbitrarily sio,i = 66, so that as 
k iterates to k = 22, ^22,2 = 5298. Table 1 shows the tabulation of the result 

5 . The number of entries at this point is only 1 3 , but a total of 3 1 entries are 
preferred, where each entry corresponds to a particular element in the pattern 

10 vector. 

The second batch of bands are obtained by taking the middle of each of the 
bands obtained in step 3. An additional 12 bands are obtained, as shows in 

Table 2: 

6. At this point there are 25 bands. The remaining six bands are obtained by 

15 combining bands from the two tables. In particular, entries 1 and 2 are merged, 

entries 3 and 4 are merged, and entries 5 and 6 are merged in both tables to 
creates six more entries, as shown in Table 3: 

Combining the above, the starting and the ending locations of the 31 bands are 
20 presented in Table 4. 



A test result on a frame of signal is shown in Figure 6 to demonstrate the robustness 
for speed changes of +/- 2% . 
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Combination of Speedup Compensation and Robust Formula 

The two methods described above for adjusting frequency band boundaries can be 
combined if speedup compensation is also incorporated. The relationship between 
5 speedup and the expansion of the frequency spectrum is exploited to combine the two 
approaches. The k-th subband, starting and the ending location = [s kl ,s k2 ] 9 has a 
robustness to speed change of +/- 2%. Each value is then multiplied by (1 + r), where 
r is the amount of speedup to [s kl , s k2 ] , followed by the roundoff method described 

above. This results in new indices [s kl9 s k2 ] whose robustness to speed change is 
1 0 shifted to r +/- 2%. Essentially, the new table is the prior Table 4, where the values 
are multiplied by (1 + 2%) and then the same roundoff method applied. Table 4 is now 
used during the registration phase to create pattern vectors from the known audio 
program that populates the database library. Table 5 is used during the detection phase 
to create the pattern vector from the detected incoming broadcast that is used in the 
15 DBS module to find matching data records in the database as described further below. 
Thus, both methods are combined. By way of example, setting r = 0.02 (2%), and 
processing every band in Table 4, a new set of subbands is calculated which is robust 
to speed change of 0 to 4%, as shown in Table 5. 

20 Table 5 is obtained with 2% speedup compensation. The new 31 pairs of starting and 
ending locations after 2% speedup compensation added to that tabulated in Table 4. 
This result is from processing the detected song from the broadcast. 

The compensation effectively positions the method to have the robustness from 0 to 
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4% speedup variations. Practitioners of ordinary skill will recognize that the same 
approach can be used to mitigate the effects of variation in the speed where the 
variation ranges above and below zero, that is, slowing down or speeding up the 
playback. 

5 

Database Search (DBS) Module 

The Database Search Module takes the pattern vector of each frame from the PG 
Module and assembles a database query in order to match that pattern vector with 
database records that have the same pattern vector. A soft matching scheme is 

10 employed to determine matches between database queries and pattern vectors stored 
in the database. In contrast, a hard matching scheme allows at most one matching 
entry for each query. The soft matching scheme allows more than one matching 
entries per query, where a match is where a pattern vector is close enough, in the 
sense of meeting an error threshold, to the query vector. The number of the matching 

15 entries can either be (i). limited to some maximum amount, or (ii) limited by the 
maximum permissible error between the query and the database entires. Either 
approach may be used. The soft matching scheme relies on the fact that the program 
patterns are being oversampled in the registration phase. For example, in the preferred 
embodiment the interframe distance used for registration is only l A of that used in the 

20 detection. Thus it is expected that if the m-th frame of a particular program is the best 
matching frame to the query, then its adjacent frames, such as (m-l)th frame and 
(m+l)th frame, will also be good matches. The combined effort of soft matching and 
sequencing schemes enhance the robustness of the detection system to varying signal 
condition inherent in the broadcasting environment. 

25 
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When matches are found, the corresponding program-id numbers and frame numbers 
in the data record is returned. The flowchart in Figure 7 illustrates the flow in DBS 
Module. Practitioners of ordinary skill in the art will recognize that a search across a 
variable to find the location of variables that match within a given tolerance in a very 
5 large database is potentially time consuming, if done in a brute force manner. In 
order to address the compute time problem, a two part search is employed. In Part 1, a 
range search scheme select those entries within a close vicinity to the query. In Part 2 
a refined search over potential candidates from Part 1 is used to select the set of 
candidates which are the closest neighbors to the query. 

10 

The steps are described in detail below: 

1. Assemble the query from the pattern vector generated by the PG Module during 
the detection phase. 

15 2. Execute a nearest neighbor search algorithm, which consists of two parts. Part 1 
exercises an approximate search methodology. In particular, a range search 
(RS) scheme is employed to determine which entries in the database falls within 
a close vicinity to the query. Part 2 exercises a fine search methodology. 
Results from Part 1 are sorted according to their distances to the query. The 

20 search algorithm can either (i) return the best M results (in terms of having 

shortest distances to the query), or (ii) return all the results with distance less 
than some prescribed threshold. Either approach may be used. As further 
described below, the nearest neighbor algorithm can be replaced with other 
algorithms that provide better compute time performance when executing the 

25 search. 
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3. If there is a match, output the program-id number and the corresponding frame 
number. If there are multiple matches, output all program-id's and 
corresponding frame numbers. 
If there is no match, output the NOMATCH flag. 

5 

Range search requires pattern vectors that match within a tolerance, not necessarily a 
perfect match in each case. From the geometrical point of view, range search 
identifies which set of the entries encompassed within a polygon where the 
10 dimensions are determined by the tolerance parameters. In the preferred embodiment, 
the polygon is a 31 dimensional hyper-cube. . 

Range Search (RS) Formulation 

In the preferred embodiment, the pattern vector is is a 1 x 31 vector: 
15 c — \c x c 2 • • • c 3 j ] , where c is the pattern vector detected where a match is sought. 
The number of bands, as described above, may be more or less than 31, with varying 
results, trading off increased accuracy for compute complexity. The search 
algorithms will be described using a 31 element vector, but practitioners of ordinary 
skill will recognize that these methods will apply with any size pattern vector. The 
20 pattern library is a M x 31 matrix, where M is the total number of pattern vectors 
stored in the database and 31 represents the number of elements in the pattern vector. 
M is a potentially huge number, as demonstrated below. Assume that the entire 
database is represented by the matrix A: 
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Those pattern vectors stored in the library are referred to as the library pattern vector. 
In the preferred embodiment, each vector z is a pattern vector of 31 elements 
5 calculated during the registration phase with known audio content for which detection 
is sought during the detection phase. During the detection phase, the identification 
exercise is to locate a set of library pattern vectors, {z_ppt}, which are being enclosed 
within the hypercube determined by the tolerance parameter. 

10 The search criteris can be represented as the identification of any z* such that 

In the preferred embodiment, LI norm is used, where ||x|| = |x 1 | + |x 2 | + — + |x 31 | is the 
LI norm ofx. Thus 




15 Here, e m , n is referred to as the nth point error between the c and z m , 

The search for z* over the entire library with the RS algorithm is based on the 
satisfaction of point error criteria. That is, each point error must be less than some 
tolerance and, in the preferred embodiment, the LI norm less than a certain amount. 
20 Practitioners of ordinary skill will recognize that the tolerance for each element and 
the LI norm may be the same or different, which changes the efficiency of searching. 



WO 2005/081829 PCT/US2005/004802 

-29- 

The determination of the tolerance is based on some statistical measure of empirically 
measured errors. Further, it is recognized that other measures of error, besides a first- 
order LI norm may be used. The search problem now becomes a range search 
problem, which is described elsewhere in the art. Reference is made to P.K. 
5 Agarwal, Range Search, in J.E. Goodman and J. O'Rourke, editors, HANDBOOK OF 
DISCRETE AND COMPUTATIONAL GEOMETRY, page 575-598, Boca Raton, NY, 
1997, CRC Press. C++ codes are also available from : Steve Skiena , The Algorithm 
Design Manual, published by Telos Pr, 1997, ISBN: 0387948600 



10 Following are the steps in the method to determine z*: 

1) Set L equal to the index set containing all the indices of library pattern 
vectors: 

Z = {l,2,3,-»,M} 

2) Start with n = 1. 

15 3) Compute e mi n between the nth element of c to the nth element of each 
Zm, n where m ranges from 1 to M. 

4) Update L to include only those indices of pattern vectors whose nth point 
error is smaller than the specified tolerance T n : 

1 < m < M, 
where 

e m9k <T k ,l<k<n 

20 T n can be set arbitrarily. In the preferred embodiment T n is set to be 10% of the 
value of c n . 

5) IfL is now an empty set AND w<31, 

Exit and issue the NO-MATCH FLAG. 
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Else: Setn = n + 1. 
Ifn>31, Go to step 6. 
Else: Go to step 3. 
6) Compute the error between all pattern vectors addressed in L to c : 
5 e m =\\z m -c\\ ; m^L 

The best solution is determined by examining all of the e m , and that will result 
with z*. Alternatively, for soft matching purposes, either of the two criteria can be 
used. Criteria 1: select only those z m with error less than some prescribed threshold 

e max* 

10 Criteria 2: select the best M candidates from L, where the M candidates are the 
least size of error to the Mth size of error. 
Once the index m with the best LI match is determined, the index is used to recover 
the data record corresponding to the pattern vector z m . The database module then 
outputs the program-id and the corresponding frame number as the output. 

15 

Note that at the start of the nth iteration, the index set L contains the indices of library 
pattern vectors whose point error from m = 1 to n-1 passes the tolerance test. At the 
start of the nth iteration, the index set L is: 



20 



\<m< M, 
L = < where 

e m , k <T t ,k = lton-l 



The flowchart of the RS algorithm is shown in Figure 8. 
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It is anticipated that the library size for application of the invention to audio 
programming, M, for 30,000 songs is in the order tens of millions. The following 
shows the calculation: 

5 Number of songs = 30,000 

Typical Song Length = 204 seconds (3 min 24 sec) 

Sampling Rate = 8,000 samples per second 

Frame Size = 1 6,3 84 samples 

Inter-Frame Distance = 4,000 samples 

10 The number of frames per song is the song length times the number of samples per 
second minus the frame size, all divided by the inter-frame distance. In the preferred 
embodiment, there are about = 404 frames 

With 30,000 songs, M = 12,1 17,120. 

With this figure, the first iteration requires around 12 million subtractions and branch 

15 statement executions to update the index set L. The next iteration will probably be 
less, but still in the order of millions. Also, memory must be allocated to hold the 
intermediate values of all of the subtraction results required for the tolerance test. 



Fast Range Search Algorithm 

20 There is an improvement to the method that minimizes the amount of subtractions that 
must be performed in order to find z*. And more importantly, the execution time does 
not scale up as fast as the size of the database, which is especially important for 
database of this size. This performance enhancement is achieved at the cost of using a 
larger amount of memory. However, practitioners of ordinary skill will recognize that 

25 because comouter memory costs have historically been reduced continuously, this is 
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r 

now a reasonable trade-off. The modification to the RS algorithm is to use indexing 
rather than computing exact error values. This modification is further explained 
below. 

5 The improved search methodology for recovering the best match between a detected 
pattern vector and pattern vectors held in the database is referred to here as the Fast 
Range Search Algorithm. As before, A is the library matrix consisting of M rows of 
pattern vectors: 
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Each row is a particular pattern vector. There are in total M pattern vectors, and in the 
preferred embodiment, each has 3 1 elements. 

15 Steps 

1. Segregate each individual column of A: 





Z h2 


Z l,31 








" Z l,2 " 




Z h31 




2 2,2 


Z 2,31 


Segregate the columns . 


Z 2,l 




Z 2,2 




Z 2,3l 


Z M,1 


Z M,2 






_ Z M,1_ 


5 




? * * *? 


_ Z M,3l_ 



20 



2. Each of the elements in the columns are sorted in an ascending order 
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Sort in Ascending order 



-2,Jt 



M,k 



f u <5 2Jt <-<f^;fc = lto31 



5 3. As a result of the sort, each element z m k is mapped to z nXk . Two cross 

indexing tables are constructed: Table Rk is a mapping of m->m 
and table T* maps m — » m , for every k= 1 to 31. 
The practitioner of ordinary skill will recognize that the sorting and table creation 
10 may occur after the registration phase but prior to the search for any matches during 
the detection phase. By having pre-sorted the pattern vectors during the registration 
phase, the system reduces the search time during the detection phase. During the 
detection phase, the method begins with a search through the sorted vectors, as 
described below. 

15 

Index Search 

Given the query vector c = [c l c 2 c 31 ] and the tolerance vector 
T = [2] T 2 ••• T 3l ] , a binary search method may be used to extract the indices of 
20 those elements that fall within the tolerance. Other search methods may be used as 
well, but the binary search, which performs in log(M) time, is preferred. 

Steps: 

1. Setk=l. 
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2. Exercise binary search to locate in the sorted column k: z A >k , m = 1 to M , the 
element iL t closest and more-than-or-equal-to c k -T k . Then exercise 
binary search again to locate the element z^ k closest and less-than-or- 
equal-to c k + T h . Thus, all the elements in the set {z^m* <m<m^} satisfy 
5 the tolerance requirement. In this manner, the binary search is used twice in 

every kth column to locate m k L and m% . 



Further, let p k be the index set containing the indices of all z A k that satisfy 
the tolerance requirement: 

10 P k -={n k L <m<m^} 



3. k = k +L if k>3 1, go to next step. 
Alternatively, the process can calculate which columns have the least number of 
bands that pass the test, and to start with that number of bands in next step. By 
15 advancing up the sorted k values where the corresponding number of bands goes from 
smallest to largest, the result can converge faster than simple increment iteration over 
k. 

4. 

Repeat steps 2 and 3 until k = 32 in order to obtain every pair of bounds: 
20 {w*,/w*},£ = lto31 , and thus determine the 31 p k 's. 

Each P k is obtained independently. For every k, all the indices enclosed 
within the pair {™*,m£},£ = lto31 can be converted back to the original 
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indices using Tk . Then, an intersection operation is run on the 31 sets of 
indices. 

An alternate way is to intersect the first two set of indices, the result is then 
intersected with the 3 rd set of indices, and so on, until the last set of indices 
have been intersected. This is the approached outlined below: 

5. Reset k = 1. 

6. Retrieve all indices in $> k and store into the array R. 

7. Use Table T k to convert all indices in R to the original indices : 



10 m — ^-^m 



Store all the indices m into a set S . 

Use Table R k+i to convert m to m : (thus the indices represented in 
column 1 are translated into their representation in column 2). Then to the 
15 results are tested to see if they are within the bound of m* +l }. 

m — 5fe±2 — >m 

Apply the tolerance test and generate 

20 R = {m 9 m k L +l <m<m^ 1 } 

In this manner, each successive p k would be the prior p k minus those 
indices that failed the tolerance test for the kth element. Thus, when k = 
30 in step 6, the p 30 are the indices that meet all 3 1 tolerance tests. 
8. k = k+l. 
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9. Go to Step 6 and loop until k = 31. 

10. Here, the set S are all the original indices after the 31 intersection 
loops. If S is empty, issue the NO-MATCH flag. Otherwise, for hard 
matching, we proceed to locate the sole winner which may be the closest 

5 candidate, for example. For soft matching, we proceed to obtain all the 

qualifying entries. 



Further speed enhancements to the fast RS algorithm 

Starting from step 4, instead of starting from k = 1, then k = 2, then k = 3, . . . , to the 
10 end, the total number of candidates in each column can be measured. The total 
number of candidates in each column is equal to the total number of candidates in 
each p k . The order of k's can then be altered so that the first k tested is where p k 
has the fewest candidates, and so on until all k's are tested. Then the order of 
intersection starts with columns with the least number of candidates. .The end result 
15 is the same as intersecting the same set of 31 indices with k incrementing 
sequentially, but by ascending the reordered k's, the number of intersecting 
operations, is reduced and thus speeds up the search. 
Search Booster: 

Practitioners of ordinary skill will recognize that the current search methodologies 
20 generally are searching on a frequency band by frequency band basis. Empirical 
studies using the preferred embodiment indicate that the initial iteration of the search 
results in 60% to 90% of the entries in the database passing the filter for that 
frequency band. Assuming a database of 6,000 song titles with 300 entries per song, 
the total number of entries to be searched is 1,800,000. With a 60% return, the system 
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has to deal with more than a million entries after the first intersection. The number of 
iterations necessary to converge on the single search result can be reduced if the size 
of the initial intersection is smaller. It is another object of the invention, referred to 
here as the booster, to pre-process the search in such a way as to reduce the number of 
5 search results in the beginning iteration of the process. 



The booster uses a different indexing scheme such that more than one frequency band 
can be lumped together. By means of the booster, a single search loop in the booster is 
equivalent to multiple loops in the range search method, and hence the search speed 
10 improved. A ranking scheme is used to determine the order of the search so as to 
minimize the number of searches for intersecting indices. To establish this rankings 
the maximum, mean and standard-deviation of the return percentile in each of the 
bands is computed during the normal range search process. These empirical results 
are used to choose which bands will be lumped together using the booster process. 

15 

The booster indexing scheme is an extension of a binary-to-decimal conversion, 
where a vector of binary-value elements is converted to a decimal integer. The 
extension is straightforward. In particular, if the base of a vector x , of size N, is M, 
where M is an integer, the conversion formula is as follows: 



20 
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£ x n if-" 

2=1 

Eqntf} 

Note that the conversion by Equation 1 is reversible, that is the equation may be used 
to convert rf- to x . Thus, the conversion possesses the one-to-one relationship so 
5 that every unique integer d 3 is calculated from a unique x . In the preferred 
embodiment, the database that houses the pattern vectors, each of the pattern element 
is stored as a 16-bit unsigned integer. This implies that each pattern vector can be 
considered as a code vector, with M = 65536 and N = 31, and a unique <i- can be 
calculated for each pattern vector. As a result of this conversion the multi- 
10 dimensional space of the pattern vectors are mapped to a one-dimensional space. 
The search for pattern vectors that are within the required distance from the query 
vector y = \y X9 y %9 ^y^ 9 referred to elsewhere as the tolerance requirement and here 
as the gap requirement, is to locate all entries x = [x 15 x 2 ,...,x w ] in the database such 
that the gap requirement \x k -y k \<Q;k = l..3l is satisfied. In the preferred 

15 embodiment, where the coding is 16 bits, the tolerance T k is 10% of the range of the 
16 bits so that Q = 10% x 64K = 6554. In practice, the value 6,000 is used. 

The booster maps the gap requirement in each band (referred to elsewhere as the 
20 tolerance requirement) to the corresponding gap requirement in d 3 . Although the 
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search can then iteratively single out all entries that satisfies all the gap requirements, 
the major difficulty of this approach is that the multiple gap requirements result in 
multiple disjoint segments on ^ . In particular, 31 iterations are required for the 
identification of the qualifying entries in where x is converted to J- , and the 
5 first loop is for band 1, the 31 st loop is for band 31. Practitioners of ordinary skill 
will recognize that by changing the number of bands in the pattern vector, the number 
of iterations would change, but the substance of the approach would be the same. 

To circumvent the technical difficulty, two compromises is made: First, only a subset 
10 of frequency bands are selected to be included in the booster, i.e., only those indices 
in the subset are coded using Equation 1. Second, a smaller base is used. The first 
compromise reduces the number of iterative loops, or specifically, the number of 
disjoint segments, so searching over every segment is practical in terms of CPU 
speed. The second compromise cuts down the memory requirement, and, more 
15 importantly, it allows for hard coding the search result of the booster (with just a 
marginal amount of RAM) to make the search within the booster very fast. 
The process for the preferred embodiment is described in detail below: 

1. SetthebaseN = 31. 

20 2. Choose 3 out of the 3 1 bands. More or fewer bands could be chosen. However if a 
large number of bands are chosen relative to the number M, then the booster method 
becomes slower and its usefulness more limited. If too few, its not accurate enough 
and does not speed up either, so an optimal number is empirically determined. In the 
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preferred embodiment, where N = 31, 3 out of the 31 are chosen. 3. This combination 
results in: 

(a) that the dynamic range of the new index is from 0 to 32767. Thus each new index 
can be coded in 2 bytes. 
5 (b) Hard-coding of the search results: Create 32768 bins: bin 0 to bin 32767. Bin m 
holds the indices of all library pattern vectors whose 3-band elements result in the 
value m after the conversion. 
4. Search Methodology: 

(a) Given a query vector j? = [y x , J> 2 >— >y n ] 
10 (b) Single out the elements in the three specified bands. 

(c) Convert the query vector using those three bands to a number using Equation 
1. 

(d) Collect all the indices of the library vectors that fulfill the gap requirement in 
the three specified bands by looking for the closest match of values m between the 

1 5 converted query and the converted library pattern vectors. 

(e) Pass the indices in (d) to the output and resume the band-by-band search 
described above on those sets of indices. 

Practitioners of ordinary skill will recognize that the conversion of the library pattern 
vectors using Equation 1 may be made prior to operation, so that the run-time 
20 computation load is reduced. 

D. Song Detection and Identification (SDI) Module. 



The SDI module takes the results of the DBS module and then provide final 
confirmation of the audio or video program identity. The SDI module contains two 
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routines: 

1 . Detection — Filtering on regularity of the detected song number: 
Irregular matches, where the DBS module returns different program-id numbers on a 
5 consecutive set of frames, is a good indication that no program is being positively 
detected. In contrast, consistent returns, where the DBS module returns consistently 
the same song number on a consecutive set of frames, indicates that a program is 
successfully detected. 

A simple algorithm based on the "majority vote rule" is used to suppress irregularity 

10 returns while detecting consistent returns. Assume that the DBS module outputs a 

particular program-id and frame-id for the ith frame of the detected program or song. Due 
to irregular returns, the result program-id will not initially be considered as a valid program 
identification in that frame. Instead, the system considers results on adjacent frames (that 
is, non-overlapping frames) of i, i+1, i+2, . . . , i+2K, where in the preferred embodiment, 

15 K is set to between 2 and 4. If there is no majority winner in these (2K + 1) frames, the 
system will issue song number = 0 to indicate null detection in the ith frame. If there is a 
winner, i.e. that at least (K + 1) frames that are contiguous to frame i produced the same 
program-id number, the system will issue for the ith frame the detected song number as 
such majority winning program-id number. Practitioners of ordinary skill will recognize 

20 that a majority vote calculation can be made in a number of ways, for example, it may be 
advantageous in certain applications to apply a stronger test, where the majority threshold 
is a value greather than K+l and less than or equal to 2K+1, where a threshold of 2K+1 
would constitute a unanimous vote. This reduces false positives at potentially the cost of 
more undetected results. For the purposes here, majority vote shall be defined to include 

25 these alternative thresholds. For computation speed, the preferred embodiment determines 
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the majority vote using a median filter. A median on an array of 2K+1 numbers, 

z = [ z i z 2 ••• z 2Jt+J ? K= 1> 2 > . . .. , is the K-th entry after Z is sorted. For example, 

if Z = [ 1, 99, 100], the median of Z is 99. The formula for such computation is stated 

below: 

5 

Assume that the DBS module returns program-id #[n] for the nth frame. To calculate the 
median for frame i: 



Let x = median([#[i] + #[i+2K$ 
1 0 Then let y = 1 - median^sgn(\#[!] - x\) sgn(j#p + 1] - x\) 
where 



sgn(|#[i + 2jq-x|)]} 



sgn(x) = 



1 x>0 
0 x-=0 
-1 jc<0 



Then, the detected result is a multiplication of x times y. . The major feature of this 
15 formula is that it can be implemented in one pass rather than an implementation 
requiring loops and a counter. 



2. Identification of programming. 



20 Given that an audio or video program is detected using majority rule, as explained 
above, the next step is to impose an additional verification test to determine if there is 
frame synchronization of the song being detected. In particular, the frame 
synchronization test checks that the frame-id number output by the DBS module for 
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each p-th frame is a monotonically increasing function over time, that is, as p 
increases. If it is not, or if the frame indices are random, the detection is declared 
void. The following are the step-by-step method of the entire SDL In cases where a 
portion of the program has been repeated, for example, in a song chorus that may be 
5 edited into the program each time, pattern vectors otherwise substantially identical but 
with varying time frames will be found by the DBS module. In these cases, the system 
carries these results along by storing them in a buffer and subjects them to the 
sequencing test explained below. As the sequencing test proceeds, some of these 
interim results will have time frame indexes that are deemed invalid under the 
10 sequencing test and will then -be ignored. Once a single interim thread survives, then 
the start and stop times of the detection are updated. 

SDI Algorithm and steps 

15 Let s p be a structure that holds the most recent 2K + 1 programid's after the p-fti 
broadcast frame has been detected: 







S p+l,l 




S p+2K,l 










S p+2K,2 


_ s pA. 




_ S +l,P 2 _ 




_ S p+2K,P 2K+l _ 



V ~lstbin ' 2nd bin (2K+l)thbin 
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Here, s mn = the n-th program id being detected in the m-th broadcast frame by the 
DBS module . Note that the P m is the size of the bin. In general, P m is different for 
different ra's. 

5 Correspondingly, f p is another structure holding the the corresponding frame 
numbers or frame indices: 



fp,\ 




fp+l,l 




fp+2K,l 


fp,2 




fp+U2 




fp+2K,2 






_fp+l,P 2 _ 






1st bin 


2nd bin 


(2K+l)th bin 



10 where f m n = the corresponding frame index of s m n . 

Also, 57 = programed of the last song or program that was successfully detected, 
such that the voting test and sequential test was successfully met. A register is created 
to hold this result until a new and different song or program is detected. 
15 Steps: 

1 . Compute the majority vote of s p 

Taking every program in the first bin of as the reference. Scan the rest of 
the 2K bins to determine if any program in the first bin pass the majority vote 
20 requirement. 
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( -\ D = Indices of entries in the first bin of that pass the majority v 

\s„„,m g D„ \ ; p 

requirement 

0 ; =0 if all the program in the first bin fail the majority vote requiren* 



If w p = 0, 

^=^+1. Go to Step 1. 
Elseif w p is a singleton (meaning a set of one element) and not equal to zero 

Set SI = w p . Go to Step 3. 
Elseif w p has more than one candidates 

Set 57 = w p (case with multiple program matches). Go to Step 3. 



10 



15 



4. 



Steps 3 to 7 are performed per s p m in w p . 

For every s m in D form a matrix A from the corresponding frame in f l 



A = 



1 

2 



f 2 



2K + 1 f %K+x 



where f t is the a frame of s p m in the /-th bin of f . 

If there is no frame in the f-th bin that belongs to s p m , f t = 0 . 

Perform the compacting of A, discarding the #-th rows in A where f q —0: 



A = 



1 

2 



/a 



2^ + 1 f 2 



2Ar+l. 



discard the qth row if f q = 0 



Cleanup A by removing rows, with the following steps: 
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A. Start with n = 1. 

B. Compute 

d t = f Ui - f ln and d 2 = k n+1 - k n . After performing step 5 by removing all 
the entries with mismatched program-id's, this step identifies only those 
5 entries that follow the sequencing correctly. 

C. Here, the quantity d x is the offset of frames between the two detected 
frames in B. This quantity can also be translated to an actual time 
offset as well: by multiplying the value by the interframe distance in 
samples and dividing by the samples per second. The quantity d 2 is 

10 the frame offset between the two broadcast frames. Now d is the ratio 

of the two offsets, representing the advance rate of the detected 
sequence. In particular, in the preferred embodiment, the system 
expects an ideal rate of 4 as the value for d. However, an elastic 
constraint on d is applied: If [d x e (4[d 2 - 1]+ 2 7 4[d 2 - 1]+ 6)], the two 

15 frames are in the right sequencing order. Thus, with d 2 - 1, an offset 

of 2 to 6 frames is expected between two adjacent broadcasting frames 
with the same program-id. If d 2 = 2, the offset is from 2+4 to 6+4 
frames. Thus the range is the same except for an additional offset of 4 
frames in the range. The values of 2 and 6 are a range centering around 

20 the ideal value 4. A range instead of a single value allows the offset to 

be a bit elastic rather than rigid. To be less elastic, one can choose the 
range to be from 3 to 5. In the same way, the range can be from 1 to 7 
to be very elastic. Go to Step D. 
Otherwise, 
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10 



15 



20 



n = n + 1, in order to sequence through all the entries in B 
Ifn<N, 

Go to Step C. 
Otherwise, 

Go to Step D. 

D. The matrix C is returned. Every row in C consists of the entries that 

satisfy the sequencing requirement. 
Compact B by deleting rows that fail to match the sequencing requirement. 
Further, note that by taking the first entry of B as the reference, if the 
second entry fails the sequencing requirement, the process can jump to the 
third entry to see if it satisfies the sequencing requirement with the first 
entry. If the second entry is satisfied with the requirement, then the second 
entry becomes the reference for third entry. 



B = 



delete rows that fail the 

-> C = 

sequencing requirement 



Jl fh 

Jp fj P 



Majority vote requirement is enforced again here. 

If the number of entries in C fails the majority vote requirement, 

the entry j is not qualified for further test, return to Step 3 for the 

next entry in D p . 

Otherwise, 

continue onto Step 6. 
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10 



20 



The majority vote test is applied again because even if the majority vote 
passes in Step 5, the majority vote test may fail after cleaning up the result 
with the sequencing rule requirement. If the revised majority vote passes, 
then a new program or song has been positively detected, otherwise, there 
is no detection. 

Let s = Number of entries (i.e. rows) in C. 

Ifs<K, 

Go to Step 9. 

Else proceed to perform regression analysis: 

A. Let C^[C n C 21 - C S J and C 2 = [C 12 C 22 - C S J be the 
first and the second columns of C respectively, where the superscript T 
denotes matrix transposition. Construct the following matrices for 
regression analysis. Regression analysis is used to calculate a linearity 
measure of the sequencing of frame-id numbers.: 







±c 2 nl 


£c nl 




2 Cnfinl 


15 




±c nl 




E = 








s 





B. Compute both the slope and the intercept 

= DE 



slope 
y - intercept 



C. Also compute the correlation coefficient r of C . 
If [r > 0.9 AND slope > 2 AND slope <6], 

the thread pertaining to the entry s p m has passed all the test and is a 
valid entry to the tracking mode. Store the entry s p m and the 
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correspondiug thread into a register called FinalList. 

Else, 

the entry s p m is discarded. 

5 Continue the test for the next entry in D p . 

8. Enter the Tracking Mode. Each thread in the Final Jist will be tracked either 
collectively or separately. 

9. Start the tracking mode: 

A. Create a small database used for the tracking: 
10 i. In the collective tracking mode, the small database contains all 

the pattern vectors of all the qualifying entries in the Final_list. 
ii. In the separate tracking mode, dedicated database containing 
just the pattern vectors for each particular entry FinaMist is 
created for that entry. 
15 B. If tracking mode = collective tracking, 

i. p=p+l: 

ii. Run detection on the G?+l)th frame of broadcast. 

iii. Update the sequence of each thread. Monitor the merit of each 
thread by observing if the thread is satisfied with the 

20 sequencing requirement. 

iv. Continue the tracking by returning to step i. if there exists at 
least one thread satisfying the sequencing requirement. 
Otherwise, exit the tracking. 



PCT/US2005/004802 



-50- 

If tracking mode = separate tracking, use dedicated database for each 
thread for the tracking. Steps are identical to that of collective tracking. 

The sequencing requirement here is the same as what is being used in 
Step 5c. That is, we expect the id of the detected frame for the new 
broadcast frame is in a monotonic increasing manner, and the 
increasing amount between successive frame of broadcast is between 2 
to 6 in the preferred embodiment. 

If for any thread being tracked, that the new broadcast failed the 
sequencing requirement relative to the previous frame, a tolerance 
policy is implemented. That is, each track can have at most Q times of 
failure, where Q = 0, 1, 2, .... If Q = 0, there is no tolerance on failing 
the sequencing requirement. 

C. After the tracking mode is terminated. Exam the merit of each thread. 
The thread that has the highest score is the winner of all in the 
FinaMist. 

i. The score can be calculated based on the error between each 
frame in the thread to the corresponding frame of the broadcast; 
or based on the duration of the thread. Or both. In our preferred 
embodiment, the duration is taken as the tracking score of each 
of thread. The one that endures the longest within the period of 
tracking is the winner thread. 

D. If multiple programs in being posted SI in Step 2. correct the posting 
by the program ed of the winning thread. 
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10. Wait for the new p-th frame from the broadcast, Go back to Step 1 . 



Practioners of ordinary skill will recognize that the values used in Step 6 for testing 
5 the linearity of the sequential frame-id's may be changed either to make the test easier 
or make the test harder to meet This controls whether the results increase false 
positives or suppress false positives while raising or lowering the number of correct 
identifications as compared to no detections. 

10 

Although the present invention has been described and illustrated in detail, it is to be 
clearly understood that the same is by way of illustration and example only, and is not 
to be taken by way of limitation. It is appreciated that various features of the 
invention which are, for clarity, described in the context of separate embodiments 
15 may also be provided in combination in a single embodiment. Conversely, various 
features of the invention which are, for brevity, described in the context of a single 
embodiment may also be provided separately or in any suitable combination. It is 
appreciated that the particular embodiment described in the Appendices is intended 
only to provide an extremely detailed disclosure of the present invention and is not 
20 intended to be limiting. It is appreciated that any of the software components of the 
present invention may, if desired, be implemented in ROM (read-only memory) form 
or stored on any kind of computer readable media, including CD-ROM, magnetic 
media, or transmitted as digital data files stored in a computer's memory. The 
software components may, generally, be implemented in hardware, if desired, using 
25 conventional techniques. 
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The spirit and scope of the present invention are to be limited only by the 
terms of the appended claims. 
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Table 1: 
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Table 3: 



603 


o43 


OA A 

o44 


1 1 C1 


1182 


1oD4 






231 / 


3243 


a A 

3244 


4541 








66 


129 




130 


255 




256 


501 




OA 

ou 


156 




157 


306 




307 


602 





Table 4: Starting and Ending Locations of the 3 1 bands in the generation of robust patterns. 
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Table 5: 
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